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Abstract 

For description of the d, s, b quark mixings the Cabibbo-Kobayashi- 
Maskawa matrices are used but they do not contain the time dependence. 
In this work the analogous matrix is obtained for the case of three neutri- 
no (z/ e , v^, v T ) mixings (oscillations) in vacuum in the general case, when 

CP violation is absent. In contrast to the quark case this matrix contains 
O . 

the time dependence. The matrix for probability of neutrino transitions 
(oscillations) in vacuum is also obtained. Naturally, it contains the time 
Oh! dependence. The matrix which does not contain the time dependence is 

obtained by using time t averaging of this matrix. Elements of this matrix 
^ ■ can be used to describe neutrino decays. 
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1 Introduction 

At present, existence of three families of leptons and quarks 
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is established [1]. In the framework of the standard model of weak inter- 
actions [2], i.e. at W, Z° boson exchanges, transitions between different 
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families of leptons or quarks with flavor number violations do not take 
place. In the quark sector, mixings between <i, s, b quarks (i. e. transitions 
between different families of quarks) are described by Cabibbo-Kobayashi- 
Maskawa matrices V q [3] 



(3) 



which are parameterized by 4 parameters, three angle 6 q , (3 q , "f q of rotation 
and one parameter of CP violation 5' . 
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c[c' 2 c' 3 - s' 2 s' 3 exp{i8' q ) c[c' 2 s 3 + s' 2 s' 3 exp(iS' q ) 
-444 - c' 2 s' 3 exp(iS' q ) -444 + c' 2 c' 3 exp(iS' q 
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4 = COs6 q , 
4 = COSf3 q , 

4 = cos7 9 , 



4 = sin6> g , 
4 = sm(3 q , 
4 = sin 7(Z , 
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Co 2 + So 2 = 1; 



(5) 



# + s 'o z = 1 



exp(^) = cos S' q + i sin ^. 

It is especially necessary to remark, that d! ', s', 6' quarks are superposi- 
tions of g?, s, 6 quarks which are eigenstates of strong interactions. 

The charged current in the standard model of weak interactions for 
three quark families has the following form [2]: 

( d\ 

3 a = (uci) Ll a V 
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and then the interaction Lagrangian is 
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From expression (7) we see that quarks are produced in the weak inter- 
actions in superposition states of d, s, b quark states - eigenstates of the 
strong interactions. And there is no time dependence there. It is necessary 



to remark that this supposition requires to be commented. This necessity 
is related with the fact that the quarks are produced in the strong inter- 
actions and then for presence of the weak interactions which violate flavor 
numbers, they are transformed into superpositions of the strong interac- 
tions eigenstates. The problem is: can the weak interactions violate aroma 
numbers at producing quarks by the strong interactions? If that's the case, 
then quarks will be produced in superposition states and there must be no 
time dependence otherwise there must be a time dependence [4]. 

The suggestion that, in analogy with K°, K° oscillations, there could 
be neutrino-antineutrino oscillations ( v — > v) was considered by Pon- 
tecorvo [5] in 1957. It was subsequently considered by Maki et al. [6] and 
Pontecorvo [7] that there could be mixings (and oscillations) of neutrinos 
of different flavors (i.e., v e — > transitions). 

Below we consider the connection between v e ,v^ v T and v\, v<i, v?, neu- 
trino states and time dependence of neutrino wave functions and transition 
probabilities at three neutrino oscillations in vacuum. 

2 Connection Between v e ,v^v T and ^1,^2,^3 Neutri- 
no States and Time Dependence of Neutrino Wave 
Functions and Transition Probabilities at Three 
Neutrino Oscillations in Vacuum 

At first we consider the connection between v e ,v^v T and ^1,^2,^3 neu- 
trino states at three neutrino transitions in vacuum and then come to 
computation of neutrino wave functions and transition probabilities. The 
case when the CP violation is absent will be discussed. 

2.1 Connection Between v e , v T and ^1,^2,^3 Neutrino States 
at Three Neutrino Transitions in Vacuum 

We can connect the wave functions of physical neutrino states ^ Ve , ^ v , ^ Vt 
with the wave functions of intermediate neutrino states ty Ul , ^ V21 \EV 3 and 
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then 
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where the neutrino mixing matrix V can be given [8, 9] in the following 
convenient form proposed by Maiani [10] (9 = 9\2, ft = #13,7 = #23): 
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c eM = cos 6, 



s eA1 = sin 0, 



c 2 +s 2 = 1- 
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(10) 
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+ S z — 1 • 



c eT = cos ft, s eT = sin ft, c 

c MT = cos7, s Mr = sin7, r 

exp(i5) = cos (5 + % sin (5, 

then Vei^iii v T neutrino states are transformed into superpositions of v\, v<i, 
neutrino states. And then the charged current in the standard model of 
weak interactions [2] for three lepton families gets the following form [2]: 



f=(enf) L l a V 
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and the interaction Lagrangian is 

C = 



9 j a W+ + h.c. 
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Using the above matrix V, we can write down it in a component-wise 
form [8, 9]: 



Vv x v^v k i 

k=l 



^ = EWn * = e,/*,T, A; = 1-3, (14) 

k=l 

where ^> Vk is a wave function of neutrino with momentum p and mass m&. 
If to separate time dependence of neutrino wave functions then we get 

* Uk (t) = e- iE ^ Uk (0), (15) 

and 

^(t) = E^%A(o)- (16) 
fc=i 

Using unitarity of matrix V or expression (14) we can rewrite expression 
(16) in the following form: 

>M*)= E E K^.e-^V^^^o), (17) 

l'=e,fj,,r k=l 

and introducing symbol b VlVll {t) 

^(^EK^e^C (18) 

fc=i 

we obtain 

E 6^W^(0), (19) 

l'=e,jjb,T 

where 6^^,(t)-is the amplitude of transition probability ^ — > 
And the corresponding transition probability ^ n — > is: 

*W*)=I E^KJ 2 . (20) 
fc=i 



It is obvious that 



E P^(t) = l. (21) 

l'=e,fi,T 

Now we will come to computation of neutrino wave functions \tV e , ^ , \I/j, t 
and a probability of transitions (oscillations) of these neutrinos. 
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2.2 Expressions for Neutrino Wave Functions and Probability 
of u e , v^, v T — > u e , Vp, v T Transitions (Oscillations) without CP 
Violation in Vacuum 

If not to take CP violation (i.e., 5 = 0) into account and use expressions 
(8-14), (14-21), then we obtain the following expressions for the amplitude 
of neutrino transitions: 

1. If primary neutrinos are v e neutrinos, then for this neutrino wave func- 
tion for v e — ► v ei v e — ► and v e — > v T transitions we get 

^e^/e^rM = [cos 2 ((3)cos 2 (9)exp(-iE 1 t)+ 

+cos 2 (f3)sin 2 (9)exp(-iE 2 t)+ 

+sin 2 {P)exp{-iE s t)}y„ e {0)+ (22) 
+ [cos((3)cos(9)exp(—iEit)(—sin( / y)sin((3)cos(9) — 

—cos(^)sin(9)) + cos((3)sin(9)exp(—iE2t)(—sin(^)sin((3)sin(9) J r 

J rcos{^)cos{9)) + sin(/3)exp(— iE 3 t)sin{^)cos{f3)]^ Uii {0) 

+ [cos(P)cos(9)exp(—iEit)(—cos('j)sin(P)cos(9) + sin( / y)sin(9))+ 

+cos((3)sin(9)exp(—iE2t)(—cos(~f)sin((3)sin(9) — sin(^)cos(9))+ 

+sin((3)exp(-iE 3 t)cos('y)cos(f3)]^ 1/T (0). 

Expression (22) can be rewritten in the following form: 

*,.- Wt (*) = KvM*vM + K^m»M + ^(*)\(o), (22') 

where 6... are coefficients before neutrino wave functions. 

1.1. Probability of v e — > v e neutrino transitions obtained from exp. (22) is 
given by the following expression: 

P„ e -„ e (*) = 1 - cos\(3)sin 2 (29)sin 2 (-t(E 1 - E 2 )/2)- 

cos 2 {9)sin 2 {2(5)sin 2 {-t{E l - E 3 )/2)- (23) 
—sin (9)sin 2 (2(3)sin 2 (-t(E 2 - E 3 )/2). 

1.2. Probability of v e — > neutrino transitions obtained from exp. (22) is 
given by the following expression: 

P Ve ^ v Jyt) = 4:Cos 2 {(3)cos{9)sin{9)[—sin{^)sin{(3)sin{9) + cos(7)cos(<9)] 
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[sin{<y)sin{P)cos{6) + cos(7)sm(#)]sm 2 (-£(£i - E 2 )/2)- (24) 

+4:Cos 2 {(3)sin{(3)cos{6)sin{^)[sin{^)sin{(3)cos{6) + cos("f)sin(9)] 

■sin 2 (-t(E 1 -E 3 )/2)-4cos 2 (P)sin(P)sin(6)sin(j)[-sin(j)sin(P)sin(6)+ 

+cos{-f)cos{6)]sin 2 {-t{E 2 - E 3 )/2). 

1.3. Probability of v e — > z/ T neutrino transitions obtained from exp. (22) is 
given by the following expression: 

P Ve ^ Vr {t) = Acos 2 {(3)cos{6)sin{6)[-cos{ 1 )sin{(3)cos{6)+ (25) 

+sin(^)sin{6)] [cos(j)sin(f3)sin(9) + sin( , j)cos(9)}sin 2 (— t{E\ — E 2 ) /2) — 

—4cos 2 (P)cos(9)sin(P)cos( / y)[—cos( / y)sin(f5)cos(6) + sin{^)sin{6)] 
•W(-t(E 1 -E 3 )/2)+4cos 2 (/?)sm(^)szn(/3)cos(7)[cos(7)sm(/?)sm(^)+ 

+sin{-f)cos{6)]sin 2 {-t{E 2 - E 3 )/2). 
The check has confirmed that P Ve ^ Ve (t) + P Ve ^ Vtl (t) + P Ve ^ Vr {t) = 1. 

2. For the case of — > z/ e , z/ T transitions we get 

^„^ e ,^ r (£) = [cos(/?)cos(0)ea;p(-i£i*) 

(—sin( r f)sin((3)cos(6) — cos{^)sin{9))+ 
+cos(P)sin(6)exp(—iE 2 t)(—sin( / y)sin(P)sin(6) + cos(7)cos(#))+ 

+sm(^)earp(-iE 3 t)sm(7)cos(^)]^ I/e (0)+ 
J r[{—sin{^)sin{(3)cos{6) — cos{^)sin{6)) 2 exp{—iEit) J r 
+(—sin( r y)sin((3)sin(6) + cos( r y)cos(9)) 2 exp(—iE 2 t) J r 
+sin 2 (j)cos 2 (f3)exp(-iE 3 t)}^^(0)+ 
[(—sin( r f)sin((3)cos(6) — cos(^)sin(6))exp(—iEit) 

(—cos( / y)sin(P)cos(6) + sm(7)sm(0))+ 
+ {—sin{^)sin{(3)sin{9) + cos( r f)cos(6))exp(—iE 2 t) (26) 
{—cos{^)sin{f3)sin{6) — sin{^j)cos{6)) + sm(7)cos 2 (/3) 
exp(-i^)cos(7)]^ r (0). 



7 



Expression (26) can be rewritten in the following form: 

**^ W M = W*)*,.(o) + W*)*^ ) + KvrWvMi (26') 

where 6... are coefficients before neutrino wave functions. 

2.1. Probability of — » neutrino transitions obtained from exp. 
(26) is given by the following expression: 

Pjr^^ft) = 1 - 4[-sm(7)sm(/?)cos(0) - cos(7)sm(6>)] 2 

[-sin(7)sm(/?)sm(^) + cos(-f)cos(9)] 2 sm 2 (-t(Ei - E 2 )/2) 
—4[—sin( / y)sin(P)cos(6) — cos( r y)sin(9)] 2 sin 2 ( r y)cos 2 (P) 

sm 2 (-t(Ei - Ed/2) (27) 
— 4[— sin( / -f)sin(P)sin(6) + cos(7)cos(#)] 2 sm 2 (7) 
cos 2 (/3)W(-*(£ 2 -£ 3 )/2). 

2.2. Probability of —> v e neutrino transitions obtained from exp. (26) is 
given by the following expression: 

Pv^vSS) = —4:[—sin( / -f)sin(P)cos(6) — cos( / -f)sin(9)}cos 2 ((3)cos(9) 

[-sin(j)sin(P)sin(6) + cos(j)cos(9)]sin(9)sin 2 (-t(E 1 - E 2 )/2) 

—4{—sin{i)sin(P)cQs(0) - cos(j)sin(9)]cos 2 ((3) (28) 
cos(9)sin(-f)sin((3)sin 2 (-t(E 1 - E 3 )/2) 
—4[—sin{ r ))sin{(3)sin{9) + cos{^)cos{9)]cos 2 {(3)sin{9)sin{^)sin{(3) 

sin 2 {-t{E 2 -E 3 )/2). 

2.3. Probability of — > v T neutrino transitions obtained from exp. (26) is 
given by the following expression: 

P v ^ Vr {t) = —4[—sin( / y)sin(P)cos(9) — cos('j)sin(9)} 

[—cos(^)sin((3)cos{9) + sin( r f)sin(9)][—sin( r f)sin(P)sin(9) 
J rCos{^i)cos{9)] [—cos("f)sin([3)sin(9) — sin( / y)cos(9)]sin 2 (—t(Ei — E 2 )/2) 
—4[—sin( r f)sin(l3)cos(9) — cos(^)sin(9)\ 
[—cos( / -f)sin(P)cos(9) + sin( / -f)sin(9)}sin( / -f)cos 2 (f5)cos( / -f) (29) 

8 



sm 2 (-*(#i - E z )/2) 
—A[—sin{^/)sin{j3)sin{0) + cos(^)cos(6)\ 
[—cos{^)sin{(3)sin{6) — sin( r y)cos(9)]sin( r y)cos 2 ({3)cos('~f) 

sin 2 (-t(E 2 -E 3 )/2). 
The check has confirmed that P v ^ Ve (t) + P v ^ v ^{t) + P v ^„ T (t) = 1. 

3. For the case of v T — > v e ,v^ v T transitions we get 

^ v T ^v„Vp, v J$) = [cos(P)cos(9)exp(-iEit) 

(—cos(^)sin((3)cos(6) + sin{^)sin(0)) + cos(/3)sin(0)exp(—iE2t) 
(—cos(j)sin(P)sin(6) — sin( / y)cos(6))+ 

+sin((3)exp(-iE 3 t)cos('y)cos(/3)]ty Ve (fy+ 
J r[(—sin( r f)sin(P)cos(6) — cos(^)sin(0))exp(—iEit) 

(—cos( , y)sin((3)cos(6) + sin(^)sin{0)) + (—sin(^)sin((3)sin(6)+ 

J rcos( r f)cos(9))exp(—iE2t)(—cos( r f)sin(P)sin(6) — sin(^)cos(9))+ 

+sin{j)cos 2 {P)exp{-iE 3 t)cos{-f)]^^{0)+ (30) 

+ [{—cos{^)sin{l3)cos{6) + sin{^)sin{9)) 2 exp{—iEit)+ 

+ {—cos{^)sin{l3)sin{6) — sin( / -f)cos(9)) 2 exp(—iE2t) + 

+cos 2 (-f)cos 2 ((3)exp(-iE 3 t)}^ l , T (0)- 
Expression (30) can be rewritten in the following form: 

^ Wr (t) = b VrVe (t)*vM + W*)*^ ) + KvM*vM> ( 30/ ) 

where 6... are coefficients before neutrino wave functions. 

3.1. Probability of v T —> v T neutrino transitions obtained from exp. (30) is 
given by the following expression: 

P„ T _>„ T (£) = 1 - 4[-cos(-f)sin(f3)cos(6) + sin{-i)sin{6)} 2 

[-cos(>y)sin(/3)sin(6) - sin{~i)cos{6)] 2 sin 2 {-t{E x - E 2 )/2) 
—4[—cos('-f)sin(f3)cos(9) + sin{^)sin{6)] 2 cos 2 {^)cos 2 {(3) 
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sin 2 (—t(Ei - E z )/2) (31) 
— 4[— cos{^)sin{(3)sin{9) — sin{^)cos{9)] 2 cos 2 {^)cos 2 {(3) 
sm 2 (-t(E 2 - E 3 )/2). 

3.2. Probability of v T — > z/ e neutrino transitions obtained from exp. (30) is 
given by the following expression: 

Pv T ^>v e (t) = —^[—cos( r f)sin((3)cos(9) + sin(^)sin(9)\ 
cos 2 {(3)cos{9)[—cos{^)sin{(3)sin{9) — sin( r y)cos(9)]sin(9) 

sm 2 (-*(^i - E 2 )/2) 
— 4[— cosily) sin(P) cos (9) + sin( / y)sin(9)]cos 2 (P)cos(9) 

cos{-f)sin{P)sin 2 {-t{Ei - E 3 )/2) (32) 
—4[—cos( / -f)sin(P)sin(9) — sin( / -f)cos(9)}cos 2 (P)sin(9)cos( / -f)sin(f3) 

sin 2 (-t(E 2 -E 3 )/2). 

3.3. Probability of v T — > neutrino transitions obtained from exp. (30) is 
given by the following expression: 

P Vr ^> v (t) = — 4[— cos( , y)sin(/3)cos(9) + sm(7)sm(0)] 

[— sin( , y)sin(f3)cos(9) — cos(j)sin(9)} 
[—cos( r f)sin((3)sin(9) — sin(^)cos(9)\ (33) 
[-sin(j)sin(f3)sin(9) + cos(7)cos(#)]sm 2 (-£(£i - £y/2) 

— 4[— cos(^)sin{(3)cos(9) + sm(7)sm(^)] 
[— sin( , y)sin((3)cos(9) — cos( / -f)sin(9)]cos( / -f)cos 2 (P)sin( / -f) 

sin\-t{Ei - £ 3 )/2) 
— 4[— cos( / y)sin((3)sin(9) — sin('j)cos(9)} 

[—sin( r f)sin(P)sin(9) + €05(7)005(6*)] 
cos(7)cos 2 (/?)sm(7)sm 2 (-£(£ 2 - £ 3 )/2). 
The check has confirmed that P Vr ^ Ve (t) + (f) + P Ur ^ UT (t) = 1. 
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We can rewrite expressions (22'), (26'), (30') for three neutrinos wave 
functions in the following compact form: 



( b VeVe {t) &„„„„(*) b VeVr {t) \ 

K^eit) 

K rVe (t) b VrVli (t) b VrVr {t) 



( <M<>) ^ 
I *, T (0) 



(34) 



Usually matrix V (8) which connects the v e , v T neutrino states with 
the ^1,^2,^3 is named Cabibbo-Kobayashi-Maskawa mixing matrix type 
for neutrinos. It is necessary to remark that in reality this matrix is matrix 
V(t) from expression (34) 

I *,„(0) ^ 
^(0) 



= V(t) 



*V(0) 



(35) 



where 



/ b VeVe (t) b VeVfl (t) b VeVr (t) \ 
V(t)= b^M W*) b^ T (t) . (36) 
V bu T u e (t) b VTV ^{t) b Vr u T (t) ) 
We see that in the case of neutrinos mixing matrix V(t) is a function of 
time t. The analogous expression (3) for quarks mixing matrix V q does not 
contain the time dependence. And at t = this mixing matrix, in contrast 
to the quark case (see (30)), has the following diagonal form: 

I 1 \ 



V(t = 0) = 



1 
\0 1) 



(37) 



We can also introduce matrix V pro b(t) probabilities of three neutrino 
transitions (oscillations) in dependence on time and write three neutrino 
transition probabilities in the following compact form: 



( Pv^vft) Pu^uM \ 

Pv^v e {t) Pv^v^t) Pv^v T (t) 



(38) 



V Pv T ^vSS) Pv^V^t) P^^t) 

If to average V(t) on time t in expression (38), then we obtain matrix 



V = V(t), 
ll 



(39) 



which has no time dependence. For this purpose it is necessary in expres- 
sions (23-25), (27-29), (31-33) for P^^t) fulfill the following replace- 
ments: 

sin 2 (—t(E 1 - E 2 )/2) = i 

sin 2 (—t(Ei - E 3 )/2) = ± (40) 

sin 2 (—t(Ei - E 3 )/2) = 1. 
Elements of matrix (39) can be used for description of neutrino decays. 

3 Conclusion 

The Cabibbo-Kobayashi-Maskawa matrices [3] are used to describe the 
d, s, b quark mixings and these matrices do not contain the time depen- 
dence. In this work the analogous matrix is obtained for the case of three 
neutrino (v e , v T ) mixings (oscillations) in vacuum in the general case, 
when CP violation is absent. In contrast to the quark case this matrix 
contains the time dependence. The matrix for probability of neutrino tran- 
sitions (oscillations) in vacuum is also obtained. Naturally, it contains the 
time dependence. The matrix which does not contain the time dependence 
is obtained by using the time t averaging of this matrix. Elements of this 
matrix can be used to describe neutrino decays. 
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